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SUMMARY 

A correlation study is made of some of the exp er mental 
data (burnout condition) with net vapor generation available m the literature. 
The data studied herein were obtained with water flow through circular tube , 
“ct^™lar cSnnels, and annular tubes. Vertically upward, horizontal, and 45 
inclined flow direct ions are considered. All channels were uniformly heated in 
the direction of flow. Correlation is achieved through modification of an ex- 
irtiS mpirtcal relation by the use of an additional channel-shape factor. The 
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INTRODUCTION 

In the design considerations of boiling fluid devices to be used for the 

?“i^!i^^ 

bolder can result in the destructive failure of the component. The heat flux 

associated with the transition from nucleate-bulk to film boiling is 

called the critical heat flux. For a flowing system this critical heat flux can 

Sa^"«cL" 

sit ion m he ~ + 1al failure The distance measured along the 

~ fom^rSf^IhrlSaftfSr— location is designated here- 
in as the critical length. 

An empirical correlation of the critical heat flux for conditions of net 
vapor generation in forced flow upward through uniformly heated circular tubes 
is presented in reference 1. Correlation was achieved for water, liquid hydr 


gen, and liquid nitrogen through the development of an empirical relation in- 
volving an enthalpy ratio and a function consisting of several dimensionless 
groups. These dimensionless groups contain variables that include mass velocity, 
tube geometry (diameter and length), and fluid property terms. The latter terms 
were evaluated at the fluid saturation temperature associated with the fluid 
pressure at the burnout location. 

In the present study, conducted at the NASA Lewis Research Center, the em- 
pirical relations for vertical upward flow through circular tubes (developed in 
ref. l) were modified and applied to data that included variations in channel 
cross section and channel inclination or orientation. The ertended correlation 
presented herein is based on sparse and incomplete ranges of data that include 
the following tube or channel characteristics: (l) vertical upward flow through 

rectangular channels and annular and circular tubes, (2) 45° inclined upward flow 
through rectangular channels and circular tubes, and (3) horizontal flow through 
annular and circular tubes. The data used for the present correlation are for 
uniformly heated tubes and channels with water as the working fluid and cover a 
saturation pressure at the burnout location that ranges from approximately 2 to 
136 atmospheres. 


DATA ANALYSIS AND PROCEDURE 

The procedures used in manipulating the data for correlation purposes con- 
sisted of a series of trfal-and-error solutions. In these solutions, the empiri- 
cal correlating relations of reference 1 were modified by the application to the 
data of successively selected dimensionless channel-shape parameters. The objec- 
tive in all such applications was to normalize as much of the noncircular channel 
critical-heat -flux data as possible on a single curve with that for circular 
tubes. Underlying this objective was the criterion that the correlation pre- 
sented in reference 1 should be considered the foundation or base upon which to 
build and should in no way be affected by the applied channel-shape parameters. 
Furthermore, any additional parameters or modified exponents to parameters should 
not yield unreasonable or unrealistic values at their limits. 

The average curve and the ±15-percent scatter-band limit curves of refer- 
ence 1 (containing about 85 percent of the data points shorn in ref. l) serve as 
the basis of comparison for the noncircular channel data. Because of the scar- 
city of data available for noncircular channels, all data points given in the 
reference reports (except as noted) are used herein, and no data selection 
similar to that employed in reference 1 was made. As in reference 1, data with 
an entering quality at the entrance to the heated channel were not included in 
the present work. 

In tab le I are listed the total number of critical-heat -flux data points 
associated with net vapor generation In the references used herein (refs. 2 to 9) 
and the range of variables covered in the references cited. It should be noted 
that for the water data considered herein the burnout location always occurred 
at the tube exit, and the pressure noted in table I Is that measured at or near 
the burnout location. 
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DEVELOPMENT OF MODIFIED CORRELATION PARAMETERS 


Channel Geometry Considerations for Vertical Upward Flow- 

In reference 1 an empirical relation was developed for circular tubes with^ 
vertical orientation relating a critical vaporization parameter C F 

eter that combined flow considerations, tube characteristics, and flu d p p 
ties. The relation developed was expressed as 


X„ = f 
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(All symbols are defined in the appendix. ) 

In applying equations (l) and (2) to noncircular channels, the diameter D 
is replaced by the equivalent diameter D e . (The equivalent diameter is defined 
as f^ tL/the channel cross-sectional area divided by the p 
heated surface. ) The use of D_ in equations (l) and (2) to put the annular 
tube and rectangular channel data on the same basis as the circular tube data 
resulted " emergence of the data. The divergence ms evidenced by values of 
S “iSer tLn the limiting value of 1.0 for the circular tubes and by a varla- 
tion of the noncircular channel data with the ratio L/D e . 

Analysis of the rectangular and annular data showed, however, that with the 
use of the defined equivalent diameter and a properly sheeted c ^me - V 

factor these data could be normalized on a single curve together with the circu 
lar tube data. A satisfactory channel-shape factor F, that depended on the 

ratio of an effective channel height to the equivalent Jf tube di- 

of an annular tube, an additional factor consisting of the ratio of the tube di 

ameters was determined by trial and error. 

The effective channel height b e for channels in which the fluid is en- 
closed or nearly enclosed by the heated perimeter (i.e., circular tubes or large 
aspect ratio rectangular channels with the two large parallel sides heated) was 
found to be a function of the radius or one-half the channel height, respectively. 
Fo^ch^nels in which the fluid is only partly enclosed by the heated perimeter 
(i e annular tubes with inner tube perimeter heated and rectangular channels 
„«h''o„“ Se side heated), the effective channel height ms found to be a func- 
tion of the entire channel height above the heated surface. 
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The channel-shape factors for circular and rectangular channels can he written as 



D 


e 


and those for annular tubes 



The channel-shape factors for the configurations considered herein are sum- 
marized in table II. 

It should be observed that as the ratio of Dp/D2 in the exponent of the 
annular-tube shape factor approaches 1.0, that is, as D-, approaches Do, the 
value of P approaches 0.5 or that of a rectangular channel. Also, for an annu 
tube having only the outer tube heated, F might be expressed as 


F = 



4 


As the exponent D^/ Dg for the aforementioned configuration approaches 1.0, F 
again approaches 0.5; however, as D^Dg approaches 0, F approaches 1.0, the 
circular tube value for F. In an annular tube configuration with both tubes 


(D /d 

heated, F might be expressed by 0.5 4 ^ , which yields values of 

range from 0.5 to 1.0 depending on the ratio D^/Dg. 


F that 


Details of the use of the channel-shape factors for data correlation are 
discussed in the following sections. 


Modification of critical vaporization parameter . - Examination of the data 
obtained with rectangular channels (ref. 2) in the high fluid quality range con- 
sistently shows values of X c exceeding 1.0 (the maximum value for circular 
tubes) by as much as 35 percent when the equivalent diameter only was used in 
equation (2). Normalization of all the data with respect to Xp (i.e. Xp not 
exceeding a value of 1.0) was obtained by the addition of an appropriate channel- 

shape factor F in equation (2) such that a more general form of X n could be 
written as ( - 


(Q /W)F°- 5 

X r = — 

C Ah 
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Modification of flow parameter GD/(l/p) * . - As in the case of the criti- 
cal vaporization parameter, the circular tube diameter D in the flow parameter 
portion of equation (l) was replaced by the equivalent diameter D e for other 
channel cross sections. For the rectangular and annular tube data, the exponent 
of the ratio L/D e , as well as the flow parameter itself, appeared to be func- 
tions of the channel-shape factor. The noncircular channel data could be corre- 
lated with circular tube data through the use of the term 


3(F) 1 ' 5 

v>J 

in place of the term (l/d) 1 ’ 5 in equation (l). For the rectangular channels 
used herein, the exponent for L/D e is a constant equal to 0.46, and the value 

of F 4 ' 5 is also a constant equal to 0.04415. (For circular tubes these values 
are 1.3 and 1.0, respectively.) For annular tubes these functions are variables 
depending on the tube-diameter ratio term in F, as discussed previously. 


The complete correlation equation for the determination of the critical heat 
flux with vertical upward flow through channels is expressed by the following 
relation: 


X C = f : 
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where X c is defined by equation (3). 

Correlation of rectangular channel data . - The correlation of critical heat- 
flux data for boiling water in terms of equations (3) and (5) for vertical upward 
flow through rectangular channels is shown in figure 1. Also shown are the 
average curve and the scatter-band limit curves obtained with the circular tube 
data of reference 1. 

The data shown in figure l(a) and some additional classified data (not 
shown herein) were used to establish the channel-shape factor and the necessary 
exponents for rectangular channels. Of the correlated data shown, 80 percent 
are within the ±15-percent scatter-band limit curves given in reference 1 for 
the circular tube data. 

The rectangular channel data shown in figure l(b) (refs. 3 and 6) constitute 
an independent check on the development of the channel-shape factor because 
(l) the data were not used in establishing the channel-shape factor, (2) the data 
were obtained from other sources, and (3) the data cover a range of larger chan- 
nel geometries than those used in establishing the channel-shape factor. Good 
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correlation of these data with the circular tube data (as well as those of the 
other channels shown (fig. 1(a)) indicates that the channel-shape factor can be 
used with a good degree of confidence to extend the correlation beyond the ranges 
specifically covered by the available data. 

Correlation of annular tube data . - The correlation of critical boiling 
heat-flux data in terms of equations (3) and (5) for vertical flow upward through 
annular tubes (refs. 5 to 7) is shown in figure 2(a), The average and scatter- 
band limit curves obtained with circular tube data are again shown for purposes 
of comparison. While the annular tube data generally lie near the circular tube 
data^ a significant and consistent trend with the amount of liquid subcooling at 
the tube entrance is exhibited in the data. In all cases the highest values of 
X c at a given flow parameter value (abscissa in fig. 2(a)) are associated with 
the largest liquid subcooling at the tube entrance, the converse also being 
generally true. No such effects due to liquid subcooling at the channel entrance 
were obtained with either circular or rectangular channel configurations. In 
terms of the parameters in equations (3) and (5), an effect of liquid subcooling 
at the tube entrance is usually associated with local-boiling critical-heat-flux 
values rather than with bulk-boiling critical-heat -flux values. On the other 
hand, the total heat input per pound of flowing liquid exceeds the enthalpy of 
the subcooled liquid, and on this basis it could be reasoned that these data 
constitute bulk-boiling critical-heat -flux values. Because the bulk temperature 
at the burnout location is not available for the data shown in figure 2(a), it 
cannot be established with complete certainty whether or not these data are 
local- or bulk-boiling critical-heat -flux values. 

An explanation for the data trend with the amount of subcooling can perhaps 
be obtained from a consideration of the two-phase flow inside annuli such as 
that used to obtain the data shown in figure 2(a). In an annulus, a progres- 
sively larger volume of cool fluid must be heated with increasing radial distance 
from the heated inner tube compared with that for a circular tube or a rectangu- 
lar channel. As a consequence of this radial increase in volume of cool fluid, 
a fluid temperature-mixing problem may exist for flow through an annular tube 
that does not exist for the other configurations. 

In an effort to correlate these same annular tube data and those for circu- 
lar tubes and rectangular channels with an analytical solution, it was found in 
reference 4 that a change in the value of the constant from 1.0 for the two 
channel geometries to 6.5 for the annular tubes was required. Reference 4 sug- 
gests that the change in the constant required in order for the equation to fit 
the data "is due primarily to retention of a relatively larger amount of the 
available liquid in the film on the unheated surfaces at the expense of the ef- 
fective liquid concentration in the channel core at the critical heat flux con- 
dition." This hypothesis, according to reference 4, is substantiated by the 
fact that annular test section data taken with a ribbed outer wall designed to 
minimize attachment of the liquid on the unheated surface are in good agreement 
with calculations" using a value of 1,0 for the constant in the equation de- 
veloped in reference 4 and by motion pictures that show "a relatively thicker, 
more placid layer of liquid flowing up the unheated wall, compared to the thin- 
ner, more agitated liquid film on the heated wall". 
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A plot of the annular tube data taken with a ribbed outer wall in the annu- 
lus (ref. 4) is shown in figures 2(b) and (c). In figure 2(b), these data are 
shown in terms of equations (3) and (5). In general, the Xq values are higher 
than those given by the average curve for the circular tube data by about 25 per- 
cent. This increase in Xq is attributed to the generally increased turbulence 
level of the two-phase fluid flowing through the annulus; turbulence has the 
effect of increasing the boiling heat-transfer coefficient, which results In the 
obtaining of higher critical-heat -flux values. Of more specific interest is the 
effect of this turbulence or better thermal mixing on the variation of Xq with 
liquid subcooling enthalpy. This variation together with that obtained by using 
a smooth outer wall In the annulus is shown in figure 2(c). For similar ranges 
of liquid subcooling enthalpy, the ribbed outer wall reduced the Xq variation 
to 12 percent from about 30 percent for the smooth outer wall. 

Further documented evidence must be obtained with annular tube configura- 
tions in order to determine whether the data shown in figure 2(a) are local- or 
bulk-boiling criti cal -heat -flux values and thereby determine whether additional 
correlation terms are required to account for liquid subcooling trends that may 
occur in some channel geometries. 


Channel Orientation Considerations 

At present, the effect of channel orientation on the critical heat flux for 
a forced flow system is difficult to evaluate because of the great lack of data. 
It can be postulated that the change in channel orientation from vertical to 
horizontal with circular tubes and rectangular channels could result in an in- 
creased agglomeration of vapor bubbles on the heated upper surface of the hori- 
zontal channel. In this event, a burnout condition would be reached at a lower 
critical heat flux for otherwise comparable flow conditions. For an annular 
tube with the inner tube heated, however, the result could be somewhat different: 
that is, the bubbles would tend to agglomerate on the upper portions of the un- 
heated tube perimeter, and thus burnout could possibly be delayed. 

Data in terms of equations (3) and (5) for a channel orientation of 45° in- 
clined upward flow (ref. 4) and horizontal flow (refs. 8 and 9) are shown in 
figures 3(a) and (b), respectively. Also shown in figure 3 are the average and 
scatter-band limit curves for vertical upward flow through circular tubes. With 
a channel inclination of 45°, the sparce data for both the circular tubes and 
rectangular channels are represented reasonably well by equations (3) and (5). 
With horizontal flow through annular tubes (fig. 3(b)), however, the data 
deviate greatly from the correlation given by equations (3) and (5). Data from 
reference 8 are generally significantly below the curves for vertical upward flow 
data, while data from reference 9 are for the most part significantly higher than 
these curves. Two preliminary data points (unpublished NASA data) with a hori- 
zontally oriented circular tube indicated reasonable agreement with the correla- 
tion given by equations (3) and (5). A data point obtained in this same study 
with a circular tube oriented vertically is also shown for comparison purposes. 

No apparent effect due to orientation for circular tubes appears evident from 
these meager data. 
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Studies of vertical upward flow through tubes and channels suggest that the 
data of reference 8 perhaps should not he considered representative of max imum 
critical-heat-flux values, hut rather of lower values associated with an unstable 
flow condition of unknown magnitude and origin with consequent submaximum 
critical-heat -flux values. Examination of the data, however, could also suggest 
that the exponents in equation ( 4 ) are a function of the orientation angle. 
Furthermore, examination of the data from reference 9 indicates that the Q c /W 
in general does not greatly exceed the liquid subcooling enthalpy at the channel 
entrance. Thus, as discussed for figure 2(a), these data are suspect in that 
they may be local-boiling critical-heat-flux values rather than values associated 
with net vapor generation. 

In view of the scarcity of critical-heat-flux data available for configura- 
tions with a horizontal orientation, consideration of an orientation function in 
equations ( 4 ) and (5) does not appear to be warranted at this time. 


CONCLUDING REMARKS 

Some of the problems encountered in correlating presently published 
critical-heat -flux data are perhaps worth emphasizing to the current researcher 
in order to obtain a greater degree of excellence in future two-phase data. All 
too often, apparent trends in the data are due to the particular research loop 
and instrumentation used by individual researchers to establish critical-heat - 
flux values. Particular attention should be exercised to avoid inadequacies and 
incompleteness in data measurements and criteria (in particular, the techniques 
and judgement used to define the often nebulous criteria of flow instability and 
its effect on two-phase heat transfer) used to establish critical-heat -flux 
values. 

It is also to be expected that critical-heat -flux values obtained in the 
future with facilities incorporating automatic burnout detectors that sense the 
large wall-temperature excursions associated with the critical heat flux may be 
higher than those reported by previous researchers. Based on limited current 
research results, these increases in critical heat flux at critical vaporization 
parameters of less than 0.7 may amount to as much as 25 percent compared with 
"manual" operation of the research equipment. The unpublished NASA data were 
obtained by using an automatic burnout detector such as mentioned above and the 
resultant values of critical vaporization parameters are considered representa- 
tive of the higher critical-heat -flux values that will be reported by research- 
ers using such devices. 


SUMMARY OF RESULT'S 

The results from a study of the critical heat flux for channels of various 
cross sections and orientation can be summarized as follows: 

1. Parameters have been evolved for boiling heat-transfer data obtained with 
water that correlate the critical boiling heat flux for forced flow upward with 
net vapor generation over a wide range of operating conditions and several 
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channel or tube configurations. The empirical relations developed herein are 
based on a previous correlation (MSA TN D-1285) that was applicable for various 
fluids, including liquid hydrogen, liquid nitrogen, and water. 

2. Limited data indicate that the correlation may also be applicable to 
inclined and horizontal channel orientations; however, in view of the scarcity 
of information available for channels with orientations other than vertical, as 
well as large variations in the results obtained by different investigators, 
further well -documented, systematically obtained data are required to verify the 
terms used in the present correlation for relating channel geometry and orienta- 
tion. 


Lewis Research Center 

Rational Aeronautics and Space Administration 
Cleveland, Ohio, April 2, 1963 
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APPENDIX - SYMBOL 


A 

b 

D 

F 

f 

G 

g 

g c 

h 

Ah 


N- 


B 


cross-sectional flow area of channel or tube, sq ft 

channel height, ft 

effective channel height, ft 

tube diameter (I.D*), ft 

equivalent diameter, 4A/ph, ft 

channel-shape factor, dimensionless 

functional notation 

mass velocity (flow rate per unit cross-sectional area of tube), 
lb mass/(hr)(sq ft) 

acceleration (gravity), 4.17X1Q 8 ft/hr 2 

conversion constant between (lb force) and (lb mass), 4.17X10 8 , 

(lb mass/lb force) (ft /hr 2 ) 

fluid enthalpy, Btu/lb mass 

total enthalpy change required for complete vaporization of incoming 
fluid, = hy. - (h^) , Btu/lb mass 


critical length from tube or channel entrance to burnout point, ft 

2 


boiling number (fluid property parameter), 


& 


P (g Or) 1 ' 5 
K v c L 


p v )g 

dimensionless 


Pr Prandtl number, dimensionless 
heated perimeter, ft 

Q total heat input from tube entrance to burnout point, Btu/hr 

q. heat flux per unit surface area, Btu/(hr)(sq ft) 

r radius, in. 

t static temperature, °P 

W weight flow rate of fluid, G(jtD /4), lb mass/hr 
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w channel width, ft 

X c critical vaporization parameter, dimensionless 

[i fluid viscosity, lb mass/ (hr) (ft) 

p fluid density, lb mass/cu ft 

a T fluid surface tension, lb force/ft 

Jj 

Subscripts: 

C critical value of heat flux (burnout condition) 

Z liquid 

s saturation 

t^ inlet temperature 

v vapor 

1 heated inside diameter (annular tube) 

2 outside diameter (annular tube) 


REFERENCES 


1- von Glahn, Uwe H. : An Empirical Correlation of Critical Boiling Heat Flux in 

Forced Flow Upward Through Uniformly Heated Tubes, NASA TN D-1285, 1962, 

2. DeBortoli, R. A., et al. : Forced-Convection Heat Transfer Burnout Studies 

for Water in Rectangular Channels and Round Tubes at Pressures Above 
500 Psia. WAPD-188, Westinghouse Electric Corp., Oct. 1958. 

5. Mendler, 0. J., Rathbun, A. S., Van Huff, N. E. , and Weiss, A.: Natural- 

Circulation Tests with Water at 800 to 2000 Psia Under Nonboiling, Local 
Boiling, and Bulk Boiling Conditions. Paper 60-HT-36, ASME, I960. 

4. Tippets, F. E. ; Critical Heat Flux and Flow Pattern Characteristics of High 

Pressure Boiling Water in Forced Convection. GEAP 3766, General Electric 
Co., Apr. 1962. 

5. Janssen, E. , and Levy, S. : Burnout Limit Curves for Boiling Water Reactors. 

APED 3892, General Electric Co.*, Apr. 14, 1962. 

6. Levy, S. : Prediction of the Critical Heat Flux in Forced Convection Flow. 

GEAP 3961, General Electric Co., June 20, 1962. 

7. Columbia University Engineering Center: Burnout Data for Internally Heated 

Annulus, Table III-5. Power Reactor Tech., vol. 4, no. 3, June 1961, p. 22. 

8. Carbon, M. W. , Gilbert, 3. E. , McNutt, C. R. , and Neider, R. : Steam Gener- 

ation in a Reactor. Problems in Nuclear Eng., vol. 1, Pergamon Press, 

1957, pp. 322-329. 

9. Foley, D* J., Batch, J. M. and McEwen, L, H. : Boiling Burnout at High Pres- 

sure in Horizontal Annuli. TID-7529, pt* 2, Reactor Heat Trans. Conf., 

1956, J. E. Yiscardi, compiler, New York (N. Y. ) Nov. 1-2, 1956. 


12 


TABLE I. - DATA SELECTION AND RANGE OF VARIABLES USED FOR CORRELATION PURPOSES 
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L Forced flow only. 

Valid data according to author. 

5 Only Includes data with fluid quality greater than 0.01. 







(a) Channel width, 1.0 Inch (Data from ref. 2). 

Figure 1. - Correlation of critical vaporization parameter for rectangular channels in terms of channel geometry 
mass velocity, and fluid property parameters. Vertical flow. 
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(b) Channel width, 1.0 and 2.1 inches. 

Figure 1. - Concluded. Correlation of critical vaporization parameter for rec- 
tangular channels in terms of channel geometry, mass velocity, and fluid 
property parameters. Vertical flow. 
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(a) Outer wall of annulus smooth. 

Figure 2. - Correlation of critical vaporization parameter for annular tubes 
in terms of tube geometry, mass velocity, and fluid property parameters. 
Vertical upward flow. 


17 



Critical vaporization parameter, 

(Q c /w)f 0 ' 5 


3 



(b) Outer wall of annulus ribbed. Tube diameters, 0.875-inch out- 
side diameter and 0.375-inch inside diameter; nominal critical 
channel length to equivalent diameter ratio, 42; fluid pres- 
sure at burnout location, 1000 pounds per square inch absolute 
(ref. 4). 

Figure 2. - Continued. Correlation of critical vaporization 
parameter for annular tubes in terms of tube geometry, mass 
velocity, and fluid property parameters. Vertical upward flow. 
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Figure 2. - Concluded. Correlation of critical vaporization parameter for annular tubes in terms 
of tube geometry, mass velocity, and fluid property parameters. Vertical upward flow. 
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(a) 45° inclined upward flow. Fluid pressure at burnout location, 2000 pounds per square 
inch absolute (ref. 2). 

Figure 3. - Correlation of critical vaporization parameter in terms of channel geometry, 
mass velocity, and fluid property parameters. Orientation other than vertical. 
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(b) Horizontal channel orientation. 

Figure 3. - Concluded. Correlation of critical vaporization parameter in terms 
of channel geometry, mass velocity, and fluid property parameters. Orienta- 
tion other than vertical. 
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